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Abstract. A high- resolution spectrum obtained with the 
low-energy transmission grating onboard the Chandra ob- 
servatory is presented and analyzed. Our analysis indi- 
cates very hot plasma with temperatures up to T ~ 15 — 
20 MK from the continuum and from ratios of hydrogen- 
like and helium-like ions of Si, Mg, and Ne. In addition 
lower temperature material is present since O vn and N VI 
are detected. Two methods for density diagnostics are ap- 
plied. The He-like triplets from N VII to Si XIII are used and 
densities around 10 11 cm -3 are found for the low temper- 
C^) . ature ions. Taking the UV radiation field from the B star 
<^> ' companion into account, we find that the low-Z ions can 
OA , be affected by the radiation field quite strongly, such that 
H I densities of 3 x 10 10 cm~ 3 are also possible, but only as- 
^ . suming that the emitting plasma is immersed in the radi- 
' ation field. For the high temperature He-like ions only low 
jy-j . density limits are found. Using ratios of Fexxi lines pro- 
{SJ ' duced at similar temperatures are sensitive to lower den- 
sities but again yield only low density limits. We thus con- 
r ~ l ' elude that the hot plasma has densities below 10 12 cnr 3 . 
, Assuming a constant pressure corona we show that the 
' characteristic loop sizes must be small compared to the 
■^j- , stellar radius and that filling factors below 0.1 are un- 
C*~) ■ likely. 

o : 

' Key words: Atomic data - Atomic processes - Tech- 
, niques: spectroscopic - Stars: individual: Algol - stars: 
' coronae - stars: late-type - stars: activity - X-rays: stars 

6 
& 



A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

missing; you have not inserted them 



ASTRONOMY 

AND 
ASTROPHYSICS 



Chandra LETGS Observation of the Active Binary Algol 

J.-U. Ness 1 , J.H.M.M. Schmitt 1 , V. Burwitz 2 , R. Mewe 3 , and P. Predehl 2 

Universitat Hamburg, Gojenbergsweg 112, D-21029 Hamburg, Germany 

2 Max-Planck-Institut fur Extraterrestrische Physik (MPE), Postfach 1603, D-85740 Garching, Germany 

3 Space Research Organization Netherlands (SRON), Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 



received February 5, 2008; accepted . . . 



1. Introduction 

Stellar coronac cannot be spatially resolved, yet they are 
thought to be highly structured just like the solar corona, 
whose X-ray emission comes almost exclusively from hot 
plasma confined in magnetic loops. So far the only way 
to infer structural information in such unresolved stellar 
point sources has been via eclipse studies in suitably cho- 
sen binary systems. Observations of the X-ray light curve 
can yield information on the location of the X-ray emis- 
sion (Prcs ct al. 1995), although the eclipse mapping re- 
construction problem is highly under-determined; after all, 
one is trying to reconstruct a three-dimensional intensity 
distribution from a one-dimensional light curve. Among a 
variety of problems discussed in detail by Schmitt (1998), 
a specific difficulty arises from the fact that in most eclips- 
ing systems both components are known or likely to be 
X-ray emitters. Obviously the reconstruction problem is 
easier to solve in those cases where one of the binary com- 
ponents is X-ray dark. At X-ray wavelengths only two such 
systems have been studied so far, the eclipsing binary sys- 
tems a CrB (cf. |Schmitt fe Kiirstci] 1993) and (3 Per (= 
Algol; |Oord fc Mewe| 1989). 

The Algol system actually consists of three compo- 
nents, a close eclipsing binary (containing a B8 main se- 
quence star and a K2IV subgiant) and a more distant 
F-type star, which is not of interest for our purposes. The 
stellar parameters of the two stars of the eclipsing binary 
system (inclination angle is i=81°) are listed in Tab. % 
Algol is one of the brightest coronal X-ray emitters in the 
soft X-ray band and has been observed with essentially all 
X-ray satellites flown so far. Particular interest in Algol's 
X-ray emission arises from the fact that no magnetic dy- 
namos and magnetic activity phenomena should occur on 
stars of spectral type B8, since such stars are fully radia- 
tive and thus the primary component of Algol should be 
X-ray dark. In consequence, all of Algol's X-ray emission 
is believed to originate from the cool secondary, which 
is rapidly rotating because it is tidally locked with the 
primary on the orbital time scale (2.8 days). We note in 



passing, however, that ther e is - in contrast to the totally 
eclipsing system a CrB (cf. [Bchmitt fc Kurster 1993) - no 
observational proof for this assumption. Nevertheless, X- 
ray eclipses at secondary optical minimum are expected, 
yet not all observations of Algol at secondary minimum 
yield evidence for such eclipses. For example, a long obser- 



vation of Algol with the EXOSAT satellite flWhite et al 



centered on secondary optical minimum showed no 
indication for any eclipse, suggesting the interpretation of 
a corona with a scale height of more than a stellar ra- 
dius or a somewhat peculiar configuration of the corona 
at the time of observation. On the other hand, a long 
ROSAT PSPC observation ( |Ottmann fc Schmitt| 1996) 
did show evidence for a partial eclipse of the quiescent 
X-ray emission, demonstrating that a significant fraction 
of the quiescent X-ray emission is emitted within a stel- 



lar radius. A BcppoSAX observation (Schmitt & Favata 



1999, |Favata fc Schmitt| 1999) of Algol showed the total 
eclipse of a long-duration flare, and a sequence of four 
ASCA observations of Algol at secondary eclipse showed 
evidence for both eclipses and absence of eclipses at dif- 
ferent occasions. 

Another method to provide information on structure 
in spatially unresolved data consists of spectroscopic 
measurements of density. If the density measurements are 
combined with the measurement of the volume emission 
measure EM, an estimate of the emitting plasma volume 
can be obtained; in an eclipsing binary these volumes 
will be subject to additional light curve constraints. With 
the high-resolution spectrometers onboard Chandra it is 
possible to carry out high-resolution X-ray spectroscopy 
for a wide range of coronal X-ray sources. We have 
obtained a Chandra high-resolution X-ray spectrum 
of Algol, which allows us to combine the information 
derived from X-ray light curves and X-ray spectroscopy. 
We will specifically discuss the Algol spectra obtained 
with the Low Energy Transmission Grating Spectrometer 
(LETGS). 



Send offprint requests to: J.-U. Ness 
Correspondence to: jness@hs.uni-hamburg.de 



J.U. Ness, J.H.M.M. Schmitt, et al.: Chandra observations of Algol 



3 



Table 1. Properties of Algol A and B: mass M, radius 
R, effective temperature T c g, logg and spectral type are 



taken from Richards (1993) and references therein. 





Algol A 


Algol B 


d/pc 


28 




M/Mq 


3.7 ±0.3 


0.81 ±0.05 


Rj Rq 


2.9 ±0.04 


3.5 ±0.1 


T cB /K 


13000 ± 500 


4500 ± 300 


logg 


4.08 


3.2 


Spectr. type 


B8V 


K2IV 



2. Instrument description and Observation 

The LETGS is a diffraction grating spectrometer covering 
the wavelength range between 5 - 175 A (0.07 - 2.5 keV) 
with a resolution A/AA ~ 2000 at the long wavelength 
end of the band pass; typical instrumental line widths 
are of the order 0.06A (FWHM) (cf., [Ness et al.| 2001a). 
A detailed description of the LETGS instrument is pre- 
sented by Predehl et al.| (1997). We note in passing that 
the LETGS uses a microchannel plate detector (HRC-S) 
placed behind the transmission grating without any sig- 
nificant intrinsic energy resolution. Thus in contrast to 
CCD based detectors the energy information for individ- 
ual counting events is solely contained in the events' spa- 
tial location. 

Accounting for the instrumental line widths the symme- 
try of the grating is sufficient to co-add both sides of the 
spectrum in order to obtain a better SN ratio ( Ness et al. 
2001a). The thus obtained spectrum is shown in Fig. |l|. 
A rich line spectrum with lines from Fe, Si, Ne, O, and 
N can be recognized between 6 and 30 A as well as many 
Fe lines above 90 A. Strong continuum emission is also ap- 
parent almost over the whole observed band pass; we note 
in passing that the spectrum shown in Fig. |l| has not been 
corrected for effective areas. 



3. X-ray fluxes and light curve 

Algol was observed with the above described instrumental 
setup between March, 12, 2000, 18:36 and March, 13, 
2000, 17:13. The total on-time was 81.41 ksec, almost 
identical to the actual exposure time. During this time 
106181 source counts were collected on the negative 
side and 101469 counts on the positive side. The high 
spectral resolution of the LETGS allows the computation 
of incident photon and energy fluxes without the need 
of any plasma emission model. In particular, since the 
LETGS wavelength range covers both the ROSAT and 
the Einstein wavelength ranges we can directly calculate 
fluxes corresponding to the respective band passes of these 
instruments without the need of any model. Using only 
bins with A c ff > 0.1cm 2 and the distance d from Tab. |] 
we compute a total X-ray luminosity of 1.4 10 31 erg/s. 



Restricting the wavelength range to the nominal ROSAT 
wavelength range (6.2-124A), we find 1.110 erg/s, 
within the Einstein band pass we find 1.0 10 31 erg/s 
(2.8-62 A). These numbers can be compared with earlier 
measurements with these instruments 
(1996) report an X-ray flux 



Bcrghofcr et al 
= 0.7 10 31 erg/s 



of L x 

measured with ROSAT. This agrees with our Chandra 
measurement to within 30% so that significant long 
term variability can be excluded. Ottmann & Schmitt 
(1996) report an X-ray luminosity of Lx — 20 10 3i erg/s 
during a flare, while their quiesc ent emission is co nsistent 
with the values reported by Bcrghofcr et al. (1996). 
Our measurement is therefore well within the range of 
luminosities found in earlier observations. 
In order to compute the ephemeris of Algol, we used 
the expression JD pnm = 2445739.003 + 2.8673285 x E 



(Kim 1989; E being an integer) for the times of primary 
minimum. Our Chandra observation covers the phases 
0.74 to 1.06, i.e., outside optical secondary minimum. In 
Fig. H we show the background-subtracted X-ray light 
curve of the LETGS data (in the ranges 10 - 120 A, 10 - 
20 A, 20 - 80 A, and 80 - 120 A). As is clear from Fig. § 
the light curve shows a more or less continuous decrease 
in intensity throughout the Chandra observations by a 
factor of 1.38 (10 - 120A), 1.22 (10 - 20A), 1.47 (20 - 
80 A), and 1.77 (80 - 120 A). Phasing of the data suggests 
the existence of a primary minimum in X-rays, when 
the late-type star is located in front of the early-type 
star, but from our discussion above this appears highly 
unlikely. Since the Algol system is known to be able to 
produce giant flares, a far more plausible interpretation 
would be to interpret the X-ray light curve as the "tail" of 
a l ong-duration flare, possibly similar to the one observed 
by [Schmitt fe Favata| (1999) with BeppoSAX. From this 
assumption, however, we would expect the radiation to 
become softer in time and to detect a cooling of the 
plasma. From Fig. ^ no evidence for softening can be 
deduced, rather the radiation becomes even harder. From 
the temperature dependent line ratios of the resonance 
lines of the H-like and He-like ions plotted in Fig. || again 
no evidence for cooling is apparent; the plasma might 
even become hotter with time, but at a very small rate. 
For oxygen the ratio raises from 3.6 ± 0.27 to 4.2 ± 
0.35, for nitrogen from 2.7 ± 0.25 to 3.3 ± 0.37 and for 
magnesium from 1.7 ± 0.12 to 2.0 ± 0.16. Thus we find 
no indication for the "tail" of a long duration flare from 
this spectral analysis. Other plausible scenarios can be 
thought of as, e.g., the time-evolution of one or more of 
the coronal active regions or rotational modulation. 



4. Data analysis 

The data extraction from the HRC-S and analysis of the 
spectra presente d in this pa per are identical to the meth- 
ods described by Ness et al. (2001a). Specifically the spec- 
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Fig. 1. Top: LETGS spectrum of Algol in the range 1-40 A. Clearly visible is the strong continuum and many prominent 
emission lines from Ne, Si, Fe, Mg, O, and N. The Ly Q lines from Mgxn (8.42A), Nex (12. 14A), OviII (18.97A) and 
Nvi (24.78A) can be recognized. 

Bottom: Wavelength range from 90-135 A. Visible are the lines from Fexvm at 94 A from Fexxi at 117.5 A and 
128.7 A, and from Fexxm/Fexx at 133A. 



tra are extracted along the spectral trace without any 
pulse height correction scheme, the background is taken 
from nearby adjacent regions on the microchannel plate. 
The two dispersion directions are co-added, but the in- 
dividual dispersed spectra can still be used to check for 



inconsistencies in the co-added spectrum. The thus ob- 
tained spectrum is shown in Fig. |l|. 
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Algol: Llghtcurve for several lines 
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Fig. 2. Light curve of Algol in a hard, medium and a soft 
energy band. 
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Fig. 3. Light curve of Algol in three different energy bands 
(top) and development of the temperature dependent ra- 
tios of Ly Q /He-like resonance line for Oxygen, Nitrogen 
and Magnesium. 



4-1- Analysis of the continuum 

The total spectrum in Fig. |l| shows two components, 
a multitude of emission lines and a significant con- 
tinuum. The shape of the continuum suggests ther- 
mal bremsstrahlung emission as the dominant contin- 
uum emission process. This assumption is supported by 
the high temperatures measured for Si, Mg, and Fe (cf. 
Tab. H). We therefore use the formula 



dW 



dV dt d\ 



2.051 x 1CT 19 nlT- 1 ' 2 



-143.9/(TA) 



A 2 



[erg cm 3 s 1 A 



(1) 

(Mewe et al. 1986) for the specific volume emissivity 
dv d dt d\ to m °del the continuum. T is the plasma tem- 
perature in K and c the speed of light in cm/s. The Gaunt 
factor G(A, T) = Gg(\,T) + G bf (A,T) is modeled as the 



sum of the free-free Gaunt factor 
fr(A,T) = 0.355A~ 06 logA 

+0.3A— log 



(2) 



0.0043 



( Mewe et alj 1986) and the bound-free Gaunt factor 
G bf (A,T) = aT»e cTd (3) 
with the p arameters a, b, c, and d adopted from 
Mewe et al. (1986). Introducing the variable x = 1/A we 
derive from Eq. [l] an expression for the total number dW of 
recorded photons divided by effective areas in each wave- 
length bin: 

dW = B * G(A, 143.9/a)/v/l43.9/a * x 2 e~ ax (4) 
with B = 2.05 x 10~ 19 EM /(ind 2 )* At* AA, a = 143.9/T, 
and x = 1/A. The plasma temperature can be derived 
from a and the emission measure EM = n 2 V from B 
while AA=0.01 A (the binsize of the spectrum), and the 
exposure time At=81.41ksec. Best fit parameters B and 
a are obtained with & x 2 n t- I n order to minimize the ef- 
fects from the emission lines, we actually use the inverse 
spectrum to be compared with \/dW{B 1 a). In each it- 
eration step we calculate the temperature from a with 
T(a) = 143. 9/ a and the Gaunt factor in each wavelength 
bin from Eqns. || and ||[ From our best fit results we find 
T = 20.7MK and EM = 67.5 x 10 52 cm- 3 . This tem- 
perature is high, we note, however, that the sensitivity to 
the temperature is quite different in the different wave- 
length ranges. In Fig. || three such wavelength ranges are 
shown for three different choices of temperature and the 
same emission measure EM = 68 x 10 52 cm -3 . Because 
of the flatness of bremsstrahlung spectra the choice of the 
temperature is not important above 18 A, but the position 
of the cutoff does depend sensitively on temperature. The 
number of emission lines between 13 and 16 A is quite high 
leading to line blends; the best fit bremsstrahlung contin- 
uum is therefore low (cf., middle panel of Fig. |4j). With 
the bottom panel we wish to show that the assumption of 
T = 14.4 MK provides a realistic description for all wave- 
length ranges and in particular for the region near the 
thermal cutoff of the X-ray spectrum. In this case 79% 
of the X-ray luminosity Lx belongs to the continuum, 
thus Lx.cont = 1-1 x 10 31 erg/sec. As a cross check we 
calculate the emission measure from the total luminosity 
and a mean radiative power loss due to bremsstrahlung of 
1.5 x 10~ 23 ergcm 3 /sec, and obtain EM = 93 x 10 52 cm -3 . 
From this exercise we note that the continuum can well 
be modeled with a bremsstrahlung spectrum, and our 
Chandra data are consistent with an emission measure 
EM = 68 x 10 52 cm" 3 and a (peak) temperature at 15 MK. 



4-2. Extracted spectra and measured line ratios 

For the analysis of emission lines we use a maximum likeli- 
hood method which compares the sum of a model and the 
instrumental background with the non-subtracted count 
spectrum. In this way Poisson statistics can be explic- 
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Fig. 4. Modeled bremsstrahlung spectrum with T — 14, 
21, and 29 MK and EM = 68 x 10 52 cm' 3 in compar- 
ison with the measured spectrum. The model is con- 
verted to LETGS counts by use of the effective areas 



from Pease ct al. (Oct. 2000). Top: very little sensitivity 
to the temperature, middle: attempt to model overlaps 
with higher temperatures, bottom: strong dependence of 
temperature suggesting 15 MK the only temperature con- 
sistent with all wavelength ranges. 




22.5 



A / A 

N VI — triplet 




A / A 

Fig. 5. (a,b) Spectrum (bold line) and best fit (thin solid 
line) for the triplets Ovu (a) and Nvi (b) for Algol. 
The dashed-dotted line represents the instrumental back- 
ground. 



itly taken into account. The model spectrum consists of 
one or more lines with a variable or fixed spacing and a 
source background, which is assumed to be constant over 
the region of interest (i.e., the spectral lines under individ- 
ual consideration). This assumption is well justified since 
both the instrumental background is quite flat as well as 
the source background (bremsstrahlung spectrum) once 
multiplied with the effective areas. The method includes 
Poisson fluctuations both in the line and background 
counts. Our code assumes the line profile functions to be 
Gaussian, but other shapes can be easily implemented. 
In this paper we used the CORA program^ version 1.2, 
for the analysis. It has been developed and described 



detailed description under 



|http : //www ,hs .uni-hamburg.de/DE/Ins/Per/Mess/Cora 



J.U. Ness, J.H.M.M. Schmitt, et al.: Chandra observations of Algol 



7 



Si xm - triplet 




A / A 
Mg XI — triplet 




X / 



Fig. 6. (a,b) Same as Fig. ||for Si XIII (a) and Mgxi (b). 
The unidentified line at 9.36 A is also fitted with 57 cts. 



by Ness et al. (2001a), and can be downloaded from 



tittp : / / ibiblio . org/pub/Linux/ science/ astronomy/ 

The analysis was performed on the basis of the 
count spectrum. The measured line counts are given in 
Tab. | and we list the best fits of the wavelengths A, 
the Gaussian line- widths, the number of line photons, 
and the source background sbg measured in counts/ A 
which is assumed constant within the individual parts 
of the spectrum under consideration. In the last column 



we list the effective areas (as provided by Pease et al 
Oct. 2000) as used for calculating line ratios needed for 
further analysis from the measurements; all errors in 
Tab. U are la errors. The first part of Tab. [| contains the 
He-like triplets Sixm, Mgxi, Neix, Ovn, and Nvi in 
combination with their H-like lines Sixiv, Mgxn, Nex, 
Ovill, and N vii. Since the Ne triplet is severely blended, 
the contaminating lines are also listed in Tab. |3J. 



For the density diagnostics of the higher temperature 
regions, five FeXXl lines are also listed in Tab. |[ together 
with the ratios of each line with respect to the FeXXl 
128.73 A line. For an estimate of optical depth effects, 
two Fe XVII lines were also measured and further analysis 
is performed in Sect. 0. 



4.2.1. Analysis of He-like line ratios 

We first discuss the lower temperature He-like line sys- 
tems from oxygen and nitrogen. In Fig. ^ (a,b) we show 
the region around the Ovn triplet at 22 A and the Nvi 
triplet at 28 A together with best fits of the resonance, 
intcrcombination, and forbidden lines. All three lines are 
clearly detected above the background, which is actually 
dominated by continuum radiation from Algol itself (cf., 
Fig. |1]). 

In Fig. | (a,b) we show the Mgxi and Sixm triplets 
together with our best fits. Obviously, the relative spec- 
tral resolution of the LETGS becomes smaller with 
smaller wavelengths, and at short wavelengths the Chan- 
dra HETGS performs far better. Still, the lines are at least 
partially resolved and line parameters can be determined 
by fitting a line template to the data whose relative posi- 
tion is fixed. In this fashion the Si xm triplet blend can be 
fitted and the determined value for the f/i ratio is consis- 
tent with the low density limit. The Mgxi triplet is more 
complicated. While the Mgxi r-line is clearly detected, 
there are no clear detections of the i and f-lincs. In par- 
ticular, the emission line feature(s) found at the expected 
position of the Mg xi f-line is unusually broad, yet we are 
not aware of other strong contaminating lines in that re- 
gion as is suggested by the fit in Fig. || (b). The determined 
line fluxes and hence line ratios do of course depend on 
the adopted background levels, yet in no case do we find 
an f/i-ratio consistent with the low density limit. Since 
this is in conflict with both the Si xm data as well as the 
Fe XXI data discussed below, we consider the "detections" 
of the Mg xi i and f lines shown in Fig. |^ (b) and reported 
in Tab. as spurious. 



4.2.2. Analysis of the Neix triplet 

The analysis of the Neix triplet is notoriously difficult, 
because of severe blending of the intercombination line 
at 13.55 A with an Fcxix line at 13.52 A (cf., Fig. ^ and 
Tab. |). The Ne ix resonance and the forbidden lines are 
clearly detected, while the intercombination line is "lost" 
in a large line blend longward of the resonance line. Our 
fits indicate 665 counts in the r-line, 385 counts in the f- 
linc, and 736 counts in the i-line blend. The question is 
how many of those 736 counts are due to the Ne IX i-line 
rather than Fexix. In the following we estimate that by, 
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Ne IX triplet with Fe XIX and Fe^jj contamination 

fit consistent with (f+i)/r=0.8 
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Fig. 7. The Ne ix triplet with blending lines from 
Fexix/xvn. The fit is constrained to (f+i)/r=0.8, i.e., 
granting a meaningfull G ratio, assuming logT < 6.7 
(5MK) for the Ne triplet. 



first, constraining the fit to G=(f+i)/r=0.8, and, second, 
by extrapolating line fluxes from other Fexix lines. 

The fit shown in Fig. [7] was performed enforcing the 
boundary condition of G=0.8. With this constraint a 
reasonable fit is obtained with 147 counts in the inter- 
combination line and the remaining 589 counts in the 
Fexix line. We now discuss whether this line count of 
589 counts is consistent wit h extrapolat ions from other 
FeXIX lines using MEKAL ( Mewe et al. 1995) for calcu- 
lating flux ratios. For this purpose we selected four other 
Fexix lines suitable for comparison, i.e., they are suffi- 
ciently isolated and/or sufficiently strong. These lines are 
located at 13.79 A, 14.67 A, 101.5 A, and 108.5 A, and our 
fit results are listed in Tab. |; for comparison we use as 
reference line the Fexix at 13.52 A. The measured flux ra- 
tios can be used for comparison with theoretical line flux 

1995). From the 
' we assume 



ratios taken from MEKAL (Mewe et al 



temperature analysis in Sect. |5| and from Tab 
a temperature of 10 MK for extrapolating flux ratios from 
the theoretical fluxes. 

The measured Fe xix and Fe xvil ratios listed in Tab. || 
indicate that the measured flux ratios are systematically 
smaller than the theoretical ratios from MEKAL. This 
effect is more significant when using the Chianti data base. 
But the general trend is quite convincing indicating the 
Fexix lines at 13.52 A and 13.79 A to be well modelled in 
Fig. ^ and in Tab. ||. Also the results for Fe xvil at 13.84 A 
seem to be realistic. 

Eventually a reduction of Fexix at 13.52 A would cure 
the difference between measured and theoretical ratios, 
leading to a higher value of the intercombination line. 
This would mean a higher G ratio (i.e., a lower temper- 



Table 2. Measured line flux ratios for Fexix and Fexvn 
lines used for consistency checks of the treatment of the 
Neix blend in Fig. 0. Values for effective areas A c g are 
taken from In-Flight Calibration by Pease et"aT| (31. Oc- 
tober 2000). The measured ratios are corrected for inter- 
stellar absorption with the values from Tab. p|. Theor etical 
flux ratios are obtained from MEKAL ( Mewe et alj 1995) 
assuming logT = 7.0 (10 MK, cf. Tab. |). 



A/A 


A [cts] 


A/A//U 




Flux ratio 






~[erg/cm 2 ] 


A/13.52 


meas. 


theor. 


Fe xix lines 




13.79 


0.24 ± 


0.30 










0.05 




13.52 


588.98 


1.66 ± 0.11 


14.67 


< 0.12 


0.11 


13.79 


147.30 


0.40 ± 0.07 


101.5 


0.07 ± 


0.16 










0.02 




14.67 


< 80. 


< 0.20 


108.5 


0.19 ± 


0.42 










0.02 




101.5 


77.670 


0.11 ± 0.02 


A/13.79 






108.5 


204.00 


0.29 ± 0.02 


13.52 


4.09 ± 


3.3 










0.85 










14.67 


< 0.49 


0.39 








101.5 


0.3 ± 


0.56 










0.08 




Fexvn lines 




108.5 


0.80 ± 


2.5 










0.18 




15.00 


1042.7 


2.56 ± 0.09 


Ratio for Fexvn A/13.84 


13.84 


114.20 


0.31 ± 0.07 


15.00 


8.16 ± 


12.9 










2.0 





ature) and a lower f/i ratio (i.e., a higher density). But 
this reduction, in combination with a reduced Fexix at 
13.79 A flux, would also require an enhancement of the 
Fexvn line at 13.84 A in order to retain the sum of the 
blend at 13.79A/13.84 A, and the ratio of this line with 
the very strong 15 A line would become smaller, which is 
the opposite effect of what would be intended with the re- 
duction. We therefore trust our results obtained from the 
constrained fit only fixing the G ratio to 0.8. 

4.2.3. Analysis of Fexxi line ratios 

The analysis of most of the Fexxi lines was straightfor- 
ward. For the Fexxi line at 121.21 A only an upper limit 
could be determined. Some difficulties were encountered 
for the Fexxi 102.22 A and the Fexxi 117.505 A lines. The 
Fexxi line at 102.22 A is partially blended with the 6 th or- 
der of Fexvn at the original wavelength at 17.054 A. The 
Fexxi line at 117.505 A is found to be very broad such 
that it is difficult to find the correct wavelength position. 
In the top panel of Fig. || our model with the isolated 
Fexix line at 101.55 A together with Fexxi at 102.22 A 
and the 6 th order of Fexvn is shown. In the bottom panel 
of Fig. the Fexxi at 117.505 A is shown in combination 
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F (101.55A) & F (102.22A) rgth d p e (17.054S) l 




101.5 



102 



102.5 



Fe 



XXI 



A / A 

7.5A 




1 17.5 



A / A 

Fig. 8. Fitting of Fexxi Top: Fexxi (102.22 A) in 
combination with Fe xix (101.55 A) and 6 th order of 
Fexvn (17.054A) Bottom: Fexxn (117. 17 A) and Fexxi 
(117.505 A) 

with the strong, isolated Fexxn line at 117. 17 A. In both 
cases the isolated lines are used to determine the line shift 
of our measurement in comparison with the theoretical 
wavelengths. In that way the expected wavelength posi- 
tion for the weaker, or blended lines under consideration 
was used for the fit. 

5. Temperature diagnostics 

We carry out temperature diagnostics using temperature 
sensitive line ratios of Ly Q /He r and of FeY/FeY + 1. 
The results are listed in Tab. || The ratios Ly a /He r 
were calculated from the line fluxes corrected for effective 
areas, as listed in Tab. |3|. We assume p lasma emiss ivities 
as calculate d in the Codes M EKAL flMcwc et al.| 1985; 
Mewe et al.| 1995) and SPEX flKaastra et al.| 1996) and 



Table 4. Temperature diagnostics from the flux ratios 
of Ly Q /Hc r for the ions Sixiv/xm, Mgxn/xi, Nex/ix, 
Oviii/vii, and Nvii/vi. T M denotes the peak formation 
temperatures. Temperatures for various Fe ions are calcu- 
lated in the same manner by using theoretical line ratios 
from MEKAL. 



Ly Q /r 



T(H-He) 
[MK] 



7m 
[MK] 



Si xiv 
Si xin 
Mg xii 
Mg xi 
Nex 
Ne ix 
O vm 
Ovu 
Nvii 
Nvi 



1.54 ± 0.14 
2.41 ± 0.31 
4.33 ± 0.26 
8.03 ± 0.71 
8.37 ± 1.28 



14.6 ± 0.5 
10.4 ± 0.5 
7.5 ± 0.2 
4.8 ± 0.2 
3.4 ± 0.2 



15.85/10.0 
10.0/6.3 
5.62/3.98 
3.16/2.2 
2.0/1.4 



flux ratio 



T/MK 



Fexvu/Fe xvni 
15.265/16.078 

Fe xxi /Fe xxn 
117.51/117.17 

Fexxi/Fexxn 
121.83/117.51 
118.66/117.51 



0.67 ± 0.07 7.93 ± 0.43 



0.19 ± 0.03 10.41 ± 0.71 



2.28 ± 0.46 
1.08 ± 0.27 



10.38 ± 0.83 
10.73 ± 1.22 



compare the measured ratios with the calculated emissiv- 
ity ratios in order to derive line formation temperatures. 
The results are listed in Tab. | as T(H-He). 

In addition we also investigate the temperature 
of Fe emitting layers with various Fe flux ratios (cf. 
Tab. H bottom). We used the photon fluxes corrected for 
effective areas from Tab. ^ and compared the ratios with 
theoretical ratios derived with MEKAL and SPEX, in the 
same manner as for the Ly a /He r ratios. The theoretical 
flux of the Fe xvni 16.078 A was corrected (enhanced) by 
a factor of 2.14 following [Mewe et al.] (2001). 
From this analysis we find a cooler component of 8MK, 
which is consistent with the Ly a /He r result for Ne. 
We also find hotter plasma at 10.5 MK in which the 
highly ionized Fe ions are formed. The ratios and derived 
temperatures are listed in Tab. ||. Clearly, a multitude of 
spectral components is present in the X-ray spectrum and 
we defer a discussion of the admissible emission measure 
distributions to a forthcoming paper. 
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Table 5. Atomic Parameters for He- like triplets. Tm is 
the peak line formation temperature (MEKAL), Ro is 
the low-density limit and N c is the density were R falls 
to half its low-density value. i?o and N c are taken from 



Pradhan fc Shull (1981) and Blumcnthal et al. (1972) 



ion 


Tm/MK 


Ro 


N C /(W 10 cm" 3 ) 


Si xiii 


10.0 


2.67 


3900 


Mgxi 


6.3 


2.6 


620 


Ne ix 


4.0 


3.5 


59.0 


Ovu 


2.2 


3.95 


3.40 


Nvi 


1.4 


6.0 


0.53 



1972, Mewe fc Schrijver 1978, Pradhan et al. 1981, 



Blumcnthal et al 



R(n e ) = - = 
i 



Ro 



He-like density from VII triplet 



1.5 

o 

o 1.0 

\ 

— 

0.5 



Pradhan fc Shull| 1981, |Pradhan| 1982, |Pradhan 1985, and 
recently Porquet fc Dubau 2000, Porquet et al. 2001, and 
Ness et al] 2001a). 

In this paper we will determine electron densities n e 
from the equation 



1 + <f>/<j> c + n e /N c y ' 

where Rq denotes the low density limit and the parameter 
N c the so-called critical density (cf. Tab. ||), around which 
the observed line ratio is density-sensitive. Finally, the pa- 
rameters 4> (the radiative absorption rate from 2 3 S to 2 3 P 
induced by an external radiation field) and cj) c describe the 
additional possible influence of the stellar radiation field 
on the depopulation of the 3 S state (cf., Sect. 3.1.2). Val- 



ues for R and N c used in this paper are listed in Tab. g. 
Also listed are the peak formation temperatures Tm for 
the ions. 



measured f/i volue 
1r/ errors 

theor. curve for 0/0 c =O 



6. Density diagnostics 

Estimates of coronal density can be obtained from the 
density sensitive f/i ratio of He-like triplets and from the 
Fexxi line ratios. The He-like Nvi, Ovu, and Neix ions 
probe the lower temperature components, while the Mgxi, 
Si xiii, and the Fexxi ions are used to probe the higher 
temperature components of the coronal plasma. The low- 
Z He-like ions are sensitive at densities log(n e ) between 
9 and 12, while the high-Z He-like ions can only be used "g 
for higher densities above log(n e ) >12. Lower densities ._ 
log(n e ) > 11 at high temperatures ~ 10 MK can be diag- ^ 
nosed from the Fexxi ratios. 



6.1. He-like ions 

6.1.1. Theory of He-like triplets 

The theory of the atomic physics of He-like 
triplets has been extensively described in the lit- 
erature (Gabriel & Jordan 1969 



0.0 k 
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1r/ errors 

theor. curve for 0/0 c =O 
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9 10 

io3(" e /cm" 3 ) 



1 1 
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Fig. 9. He-like densities n e from the ions Ovu, Nvi, 
and Neix. For comparison the theoretical curve assum- 
ing (f>/<j) c = is also plotted (cf. Eq. ||). 
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measured f/i value 3.66 
Iff errors 



0/^=0.000 
N C =3.900E+13 
R„=2.670 
f/l=3.66 



42.90 



8 10 12 
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Table 6. Investigation of the influence of the stellar radi- 
ation field originating from the B8 star. Measured fluxes 
from the IUE satellite Fx are converted to intensity I\ 
taking into account limb darkening effects using e — 0.44 



for the U band (Di'az-Cordovcs ct al. 1995) 



Nvi 



Ovu 



Ne ix 



Am/A 
WeV 

Ix 



1900 
552.1 
20 ±2 

9.1 ±0.91 



1630 
739.3 
22 ± 2.2 

9.98 ±1.0 



1266 
1195.3 
20 ±6 

1 ±2.7 



(1 7 ^iip ) 

v cm A 5 A utrd ' 

T rad /K 12066 ± 186 12708 ± 190 13686 ± 560 
dilution factor 0.01 

24.74 ± 2.41 2.18 ± 0.29 0.1 ± 0.03 



Fig. 10. He-like densities n e from the Sixni ion. 



6.1.2. Influence of the stellar radiation field 

Given the effective temperature of 13000 K of Algol A and 
its close proximity to Algol B (cf., Tab. 0), we must check 
to what extent the X-ray radiation originating from the 
corona of Algol B is influenced by the UV radiation from 
Algol A; the UV-radiation from the Algol B component it- 
self is small when comparing its effective temperature with 
Capella (T e g <5000K) and the values computed for cj)/(f> c 
to be used in Eq. |^ as derived by Ness et al] (2001a). Also 
the measurements of other G and K type stars, as pre- 
sented by [Ness et aT| (2001c) suggest Algol B not to con- 



tribute to the total radiation field. Following Ness et al 



(2001a) we used IUE measurements of Algol in order to 
derive radiation temperatures for the desired wavelengths, 
listed in Tab. [| From this we calculated values </>/0 c for 
N vi, O vii, and Ne ix using a dilution factor of 



W 



1/2- 



0.01 



(6) 



(Mewe & Schrijver 1978) using as distance between 
Algol A and Ba = 14.6 Rq ( |Forbes| 1997) and the radius 
of Algol A r* = 3.5 i?« (Tab. |l|). The result of our 
analysis is listed in Tab. 0. Since we assumed the worst 
case scenario, we find the maximally possible values <p/^) c , 
so that the effects are not negligible for O vii and not 
zero even for Neix. We do point out, however, that in 
the phases between 0.74 and 1.06 (cf. Sect. ||) much of 
the visible coronal emission could originate from regions 
not illuminated by the B star. With an inclination angle 
of 81.5° and the assumption of a uniform distribution of 
the coronal plasma, geometrical considerations lead to 
the result that at phase <f> = 1 only plasma near the polar 
regions can be illuminated by the primary component. A 
detailed discussion of these geometrical considerations is 
given in Ness et al. (2001b). In the following analysis we 



calculate coronal densities with and without the effects 
of the radiation field from the B star, i.e., assuming 
= and with the values for 4>/4> c from Tab. o. 



6.1.3. Densities with the He-like ions 



The measured ratios f/i, corrected for A e g (Pease et al 
Oct. 2000; the values are listed in the last column of 
Tab. |), as quoted in Tab. |, were used for density diag- 
nostics. The theoretical curves from Eq. [|, with the values 
4>/4> c from Tab. |^ and the other parameters from Tab. ^, 
are plotted in Figs. ^ and |l^ for each ion in compari- 
son with our measurements for f/i with ler errors. For the 
low-Z elements O, N, and Ne we also considered the case 
of no radiation field from the B star affecting the emit- 
ting layers with a line-dotted line in Fig. ||. The influence 
is most severe for Nvi and Ovu, but is also visible for 
Ne ix. Assuming a negligible radiation field from the pri- 
mary component we obtain definite deviations from the 
low density limit for N vi and O vii and a marginal devi- 
ation for NeiX; in the latter the latter case the 2a error 
includes the low density case. In these cases we find densi- 
ties between 1-2 xl0 11 cm -3 (cf., Tab. |?]). Assuming the 
full radiation field supplied from the B star to be effective 
the sensitivity of our measurements to detect densities is 
significantly reduced because the difference between the 
R-value in the high density case (R = 0) and low-density 
case becomes smaller and smaller. Specifically, for our Al- 
gol LETGS spectrum we find that the data are consistent 
with the low-density limit for nitrogen and neon, and even 
for oxygen the low density limit is included within the 2a 
error bars. For silicon and magnesium radiation effects are 
unimportant. The derived f/i-ratio for silicon is consistent 
with the low density limit, while our measurements for 
magnesium (forma lly) y ield densities of w logn e = 13; as 
discussed in Sect. 4.2.1 and Fig. || (b), we consider the 
measurements of the i and f lines in magnesium spurious. 
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6.2. Density diagnostic with Fexxi 

In addition to He-like ions, ions with more than two elec- 
trons can be used as a density diagnostic. The ground con- 
figuration of Fexxi ls 2 2s 2 2p 2 splits up into 3 P, 1 D, 1 S. 
The energy difference between the ground state 3 Po and 
the excited levels 3 Pi and 3 P 2 is 9 eV and 14.5 eV, respec- 
tively, and 30 eV between ground state and 1 D 2 . In low- 
density plasmas virtually all atoms are in the ground state 
3 P, while in high-density plasmas a Boltzmann equilib- 
rium with the higher level will be obtained. Consequently, 
from excited levels certain lines will appear only in high- 
density plasmas. In contrast to He-like lines the appear- 
ance of certain lines is an indicator of high-density plas- 
mas. 

Our measurements of four Fexxi ratios, corrected for 
effective areas and interstellar absorption are listed in 
Tab. H (effective areas and values from interstellar ab- 
sorption are listed in the last two columns) and plotted 
in comparison with t heoretical fl ux ratio vs. n e curves ob- 
tained with Chianti ( Dere et al. 2001) in Figs, [ll] and |l2|. 
As can be seen from Figs. |ll|and |l2] , only low density lim- 
its or upper limits are obtained for all Fexxi line ratios. 
Similar results are obtained when using the ratios from 
Brickhouse et al. (1995). The most sensitive upper limit 
(logn e < 11.52) comes from the Fexxi 102.2 A/128.7 A 
ratio, and is a factor ~ 24 below the upper limit derived for 
the Si xin triplet, which is formed at similar temperatures. 
This comparison clearly shows that Fe xxi line ratios yield 
far more sensitive density constraints at high temperature 
as compared to He-like triplets from magnesium, silicon, 
and higher ions. 



1 .00 



Algol: FE XXI 97.87/128.73 



7. Optical depth effects 

The Chandra Algol spectrum contains a number of strong 
Fe xvii lines which are sensitive to optical depth effects. In 
order to estimate the optical depth, we use an "escape fac- 
tor" model with a homogeneous mixture of emitters and 



absorbers in a slab geometry (e.g., Kaastra & Mewe 1995, 
Mcwc et al.| 2001). In this geometry the escape factor is 
P(t) ~ [1 + 0.43-r] -1 . Significant optical depths will lead 
to resonant scattering in strong lines. 

If one considers two lines produced by the same ion, 
the ratio of the optical depths is given by 

- = FT- ^ 

Therefore, one must look for lines from the same ion and 
large difference in oscillator strength or wavelength. In the 
solar context one usually studies the ratio of the strong 
Fexvn 15.03 A resonance line (/ = 2.66) and another, 
nearby Fexvn line with a small oscillator strength (e.g., 
15.265 A, / = 0.593). Estimating the escape factor P from 



P : 



15.03 A 
\ 15.27 A 

tical depth. 



/ 



15.03 A 
15.27 A 



we can deduce the op- 
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Fig. 11. Density diagnostic with Fexxi ions, lcr errors 
are indicated with dashed line style (used 2<r error for 
102 A ratio). 



For the measured Fexvn 15.03/15.265 photon flux ra- 
tio we obtain a formal fit result of 2.81 ± 0.25 (cf. Tab. ||). 
This ratio does depend on the assumed background value 
and can vary between 2.32 and 3.45 when varying the 
source background (between 4000 and 5000 cts/A; comp. 
Tab. |^). In this particular wavelength region (cf., Fig. |[ 
middle) line blending is severe and a correct "eye" place- 
ment of the continuum is difficult. However, our contin- 
uum modeling predicts relative stable values of 4500 cts/A 
so that we are confident that our quoted result is correct 
and not affected by systematic errors from the placement 
of the continuum background. We point out that this mea- 
surement agrees remarkably well with the same line ratio 



as measured in Capella (Mewe et al. 2001); this is interest 



ing because in Capella no measurements of density could 
be obtained. 

The measured photon flux ratio must be compared to 
the r = flux ratio, which can be deduced either from 
theory or laboratory measurements. With SPEX we pre- 
dict a ratio of 3.5, thus P = (2.81 ± 0.25) /3.5 = 0.8±0.07 , 
and therefore r = 0.57 ± 0.26. With Chianti (Dere et al 
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Fig. 12. Density diagnostic with Fexxi ions. 



2001), a ratio of m 4 is expected for r = 0, and using 
P = (2.81 ± 0.25)/4 = 0.7 ± 0.06 we find r = 0.98 ± 0.3. 
In either case we find optical depths significantly different 
from zero (at ^ 2 — 3cr level) . 

Unfortunately theory does not agree with experiment. 
The very same line ratio can be measured in the Liver- 



more Electron Beam Ion Trap (EBIT; Brown ct al. 2001, 
Laming et al. 2000). These experiments typically yield 
Fexvn 15.03/15.265 photon flux ratios in the range 2.5 
-3.0, which are significantly different from those expected 



theoretically. Also, Brown et al. (2001) point out that con- 
tamination of the 15.265 A with Fexvi further lowers the 
observed 15.03/15.265 photon flux ratio. Comparing the 
Algol (and Ca pella) flux rat ios in 15.03/15.265 to the val- 
ues quoted by Brown et al. (2001) we therefore conclude 
that the observations are fully consistent with an opti- 
cal thin plasma without any significant optical depth with 
possibly some contamination arising from Fe xvi. It is wor- 
rying that the theoretically predicted emission from some 
of the strongest emission lines observed in solar and stellar 
X-ray spectra appears to be wrong by « 30 percent. 



Table 7. Summary of derived densities. 

* Resulting density values accounting for stellar radiation 

fields. 

** 2a upper limit. 



He-like triplets 





f/i 


log(n e ) 


Si xiii 


3.66 ± 1.44 


<12.9 


Mg xi 


0.90 ± 0.37 


n.a. 


Neix 


2.62 ± 0.61 


11.30 ± 0.60 


*Neix 




<11.57 


Ovn 


0.94 ± 0.20 


11.04 ± 0.13 


*Ovn 




10.54 ± 0.53 


Nvi 


0.21 ± 0.09 


11.16 ± 0.23 


*Nvi 




10.13 ± 0.93 


Fexxi ratios 


A/A 


ratio 


log(n e ) 


97.87 


0.12 ± 0.03 


<12.19 


102.22 


0.18 ± 0.03 


< 11.52** 


117.51 


0.20 ± 0.04 


<12.63 


121.22 


<0.09 


<12.4 



8. Discussion and Interpretation 

From the derived parameters we can obtain structural 
information on Algol's corona. From the lower temper- 
ature He-like ions definite coronal densities could be 
determined, from the higher temperature He-like ions 
as well as from the Fexxi line ratios (cf., Tab. 0) upper 
limits to the coronal density at a temperature of 10 MK 



can be derived. Since EM 



V, typical coronal volumes 



V can be estimated from the recorded emission measures. 

We now assume that Algol's corona is composed of a 
multitude of individual but identical loops, all of which 
obey the loop scaling equation (Rosner et al. 1978) 



,L = 1.3 x 10 T, 



apex ' 



(8) 

with the densities n e , the apex temperature T apex , and the 
loop semilength L, all in cgs units. Of course, the apex 
temperature T apox cannot be directly measured, it must 
however be close to the formation temperature of the ion 
with the highest ionization temperature and it can cer- 
tainly not exceed the temperature of the bremsstrahlung 
continuum (Sect. |4.l| ). It is also clear that in Algol one is 
most likely dealing with a distribution of X-ray emitting 
loops as is the case in the solar corona. However, our data 
docs not allow us to constrain the parameters of these 
distributions and therefore we model the distributions by 
their presumed mean values. From Tab. |] we conclude 
that the apex temperatures T apex must be approximately 
15 MK, which is also supported by the continuum temper- 
ature, while from Tab. |we conclude that logn e < 11.5. 
Obviously the densities at 15 MK could be arbitrarily low. 
However, if we assume that the material at 15 MK is in 
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thermal pressure equilibrium with the material at 2 MK 
(where the Ovn line is produced) we estimate densities 
in the range logn e ~ 9.7 - 10.2, depending on whether 
the X-ray emitting corona of Algol B is immersed in the 
UV radiation field of Algol A or not. For logn e = 11.5 
we compute L = 9 x 10 8 cm while for densities logn e = 
9.7 - 10.2 we find L = 180 - 460 x 10 s cm. In either case, 
the loop lengths are much smaller than the stellar radius 
of 2.45 xl0 11 cm so that we assume to deal with an es- 
sentially planar geometry. If we assume "canonical loops" 
with circular cross sections of a tenth of the loop length, 
the volumes of such coronal building blocks are in the 
range 2.510 25 cm 3 to 3 1 30 cm 3 , again depending on the 
assumed density, and in all cases one requires at least 1000 
loops (or more) in order to account for the total observed 
emission measure. 

Let us next assume that the considered loops are semi- 
circular and extend to height 

E-^. (9) 
We define the coronal filling factor as the ratio between 
the available volume 4ttR 2 H and the actual volume of X- 
ray emitting material V COI = EM/n 2 : 



f = 



V a 



EM 



(10) 



^available 8n 2 i? 2 L ' 

where R denotes the stellar radius. Inserting the scaling 
law to replace L with T ap ox and n e , we finally find 

/ = 9.6x10-^^ = 9.6-^ (ii) 



n R 2 T 2 

"e 11 ^apcx 



ni R 2 n Ti 

[EM 52 = EM/W 52 cm-'^Rn = i?/10 n cm, 

T 6 = T apcx /10 6 K,mo = n e /10 10 cm" 3 ]. 
Using the continuum values EM52 = 68 and Tq = 15 we 
find / x n e — 5 x 10 9 cm -3 . Since the filling factor / can be 
unity at most, we deduce that n m in — 5 x 10 9 cm -3 and 
the actual density ought to exceed that value. Interest- 
ingly, n m i n T apex /Tovn is close to the density determined 
from demanding that the high temperature material is at 
the same pressure as the O VII emitting material (assum- 
ing that it is immersed in the primary's radiation field). 
Alternatively, density values of 1.5 x 10 10 cm -3 lead to fill- 
ing factors of « 0.3. These densities can only be reached 
without the illuminating UV radiation field of the pri- 
mary component. The overall phase range covered by the 
Chandra LETGS observations is 0.32. During that period 
the star rotates by 115°; therefore most of the hemisphere 
visible at the beginning of our observations are not visible 
at the end. The observed light curve (cf. Fig. |) with its 
slowly decreasing trend does not appear to give the im- 
pression of coming from a number of unrelated regions, 
rather, it appears to come from the same region on the 
star. This is only possible if the region is located near 
the pole (in which case it would be immersed in the radi- 
ation field) or if it is fortuitously located in a equatorial 
region which happens to be visible during the whole phase 
interval. While the latter possibility cannot be excluded, 
the former option appears more plausible given the previ- 



ous evidence for polar activity in Algol ( Schmitt & Favata 
1999), the available VLBI images of Algol QMutel et al 



1998) and the general occurrence of polar spots in Dopplcr 
images in general. It is attractive to interpret the observed 
X-ray emission as arising from the circumpolar regions, 
possibly from a long-duration flare similar to the one ob- 



served by Schmitt & Favata (1999). Assuming as height 
the limit inferred from the eclipsing BeppoSAX flare (H 
< 7 10 10 cm) and as filling factor arbitrarily a value of 1 
percent (which is the surface area of the circumpolar polar 
cap), we find a volume of 4 x 10 32 cm 3 , which requires den- 
sities of 2.5 x 10 10 cm -3 when combined with the observed 
emission measures. These densities refer to the hot mate- 
rial at temperatures of « 15 MK. Clearly, these densities 
are fully consistent with the observed upper limits. 
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Fig. 13. Temperatures required for asserting certain filling 
factors (cf. Eq. [l^) assuming constant pressures calculated 
from the densities and temperatures measured for Ovn. 
For the emission measure we use the value obtained from 
the continuum EM = 68 x 10 52 cm~ 3 . 



Continuing in our picture of "canonical" loops, we 
can compare a typical pressure scale height H p = 
5OOOT apox /( 5 / 50 ) = 3.8 x 10 n cm using T apcx = 15 MK 
and g/gQ = 0.2, with a typical loop size L « 30 x 10 s cm 
(Eq. ||). Since H p >> L, the assumption of constant pres- 
sure is well justified. Hence P — n e * T must be constant 
and thus Eq. [ll] can be rewritten as 

/ = 9 - 6 /f% <=► T 6 = 9.6^- (12) 

[Pie =n e *T/10 16 cm- 3 K], 
which relates the unknown filling factor with the - in prin- 
ciple - measured parameters stellar radius, emission mea- 
sure and coronal pressure. In Fig. |l3| we plot the depen- 
dence of temperatures on filling factor for different pres- 
sures, using Eq. [L2t the emission measure obtained from 



the continuum (Sect. 4.1) was used. The pressures are 
calculated from the measured temperatures (Tab. ||) and 
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densities (Tab. [?]) for Ovn. In the case of the radiation 
field from the B star to influence the density diagnostics, 
a lower limit for the pressure of P = 8 x 10 16 cm~ 3 K, 
using the peak formation temperature T = 2.2 MK is cal- 
culated, while we obtain P = 50 x 10 16 cm~ 3 K for the 
higher density from omitting the possible effects from the 
radiation field and the higher temperature obtained from 
the Ly Q /r ratio (Tab. g). 

It is clear that T must exceed 10 MK and is very 
likely below 30 MK. For low pressure (Pi6 = 8) the filling 
factor would have to exceed 0.45, while for high pressure 
(Pi6 = 50) the filling factor would have to be in the 
interval 0.07 < / < 0.21. Since we feel that filling factors 
/ ~ 1 are unlikely, we favor a high pressure scenario. 



9. Conclusions 

The high-resolution spectrum obtained with the LETGS 
on board the Chandra observatory shows a large number 
of emission lines and an unusually strong continuum. 
We analyze both the emission lines and the continuum 
relying exclusively on ratios of individual lines. The 
LETGS observation interval covered the orbital phases 
between 0.74 to 1.06, the light curve appears to indicate 
a relative state of quiescence. A slow decline is seen, but 
cannot be attributed to the decay of a giant flare, since 
no cooling and no softening of the overall emission is 
detected. The total luminosity is Lx = 1.4 10 31 erg/sec, 
roughly consistent with the X-ray luminosities found 
previously with Einstein and ROSAT. 

We analyzed the continuum in order to determine an 
upper temperature and an overall emission measure. The 
continuum can well be modeled with a bremsstrahlung 
continuum and a temperature of 15 MK can be derived 
which is consistent with the peak formation temperature 
of, e.g., Si xiv which is clearly detected. Optical depth 
effects leading to resonant line scattering were analyzed 
and ruled out. The ratio of Fexvn lines at 15 A and 
15.27A are used for testing optical depth effects, but 
inconsistencies in the atomic line data were found. From 
our findings in comparison with other measurements with 
the LETGS we conclude that a line ratio of ~ 2.8 must 
occur in plasmas with r = 0. 

The measurement of plasma temperatures is carried 
out using line ratios of ions from the same element 
in adjacent ionization stages as, e.g., Ly Q and He-like 
resonance lines. At least two temperature components 
are found, allowing Nvi (3.2 MK) and Sixiv (14.2 MK) 
to be formed. This is consistent with our findings for 
Fexxi/Fexxn (« 10 MK). 

For the density diagnostics with the He- like triplets the 
UV radiation field originating from the B star companion 



was analyzed. We found significant effects for the ions 
Nvi, Ovn, and even for Neix, however, the illumination 
geometry of the primary B-type star is unclear. We there- 
fore considered both the case with full illumination and 
no illumination of the coronal plasma. The determined 
densities range around 10 10 — 10 n cm~ 3 , while for the 
ions Si XIII and NelXwe find only upper limits. This is 
also found for all Fexxi line ratios, but these upper limits 
are consistent with densities 10 10 — 10 11 cm -3 as well. 

A detailed analysis discussing structural information 
deduced from the derived temperatures and densities 
is presented. The assumptions used for this discussion 
are that Algol's corona is composed of a multitude of 
individual but identical (semi circular canonical) loops, 
all of which obey the RTV loop scaling equation. We 
further assume that the material is in thermal pressure 
equilibrium at 15 MK. From these assumptions we 
conclude from our measurements that these loops have 
semi- lengths of L = 9 x 10 8 cm. The assumption of the 
emitting plasma to be immersed in the external radiation 
field leads us to lower densities and thus loop lengths of 
up to L = 460 x 10 8 cm, which is, however, still much 
smaller than the stellar radius. The emission measure 
derived from the analysis of the continuum in combina- 
tion with the densities allows us to define constraints on 
the coronal filling factor. However, different pressures 
ranging between 8 and 50 x 10 16 cm~ 3 K are possible, 
depending on whether the emitting plasma is immersed in 
the external radiation field or not. A low pressure plasma 
(i.e., effects from the radiation field are severe) requires a 
large filling factor, while the high pressure scenario allows 
a lower filling factor. 
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Table 3. Measured line counts for Algol with ler errors. Values of the total effective areas A e g are taken from In- 
Flight Calibration by [Pease et al.| (31. October 2000). The Transmissions from the ISM are based onN(Hi) = 2.5xl0 18 , 
N(Hei)/N(Hi) = 0.09, N(Hen)/N(Hi) = 0.01 





A [A] 


a [A] 


A [cts] 


sbg 
[cts/ A] 


Robs —f/i 

[R]& 


<J"obs — — 

[Gp 


A cff 
[cm 2 ] 


He-like and H-like (cf. Sects. 


^ and 


6.1.,i) 








Si xiv 
Sixnir 

f 


6.19 ± 0.0016 
6.65 ± 0.004 
6.69 ± 0.014 
6.74 ± 0.006 


0.025 ± 0.002 
0.022± 0.003 


658.32 ± 35.06 
480.7 ± 35.95 
86.4 ± 33.10 
314.1 ± 30.60 


5119 
5300 


3.64 ± 1.44 
[3.66 ± 1.45] (2) 


0.83 ± 0.15 
[0.84 ± 0.15] (2) 


36.22 
37.54 
37.38 
37.16 


Mg xii 
Mgxir 

f 


8.42 ± 0.0016 
9.17 ± 0.0015 
9.23 ± 0.0015 
9.31 ± 0.0015 


0.02 ± 0.002 
0.02 ± 0.002 


578.04 ± 33.13 
224.12 ± 26.15 
84.85 ± 23.29 
76.09 ± 22.80 


5900 
5400 


0.90 ± 0.36 
[0.90 ± 0.37] (2) 


0.72 ± 0.22 
[0.73 ± 0.23] (2) 


32.28 
27.69 
27.42 
27.24 


Nex 
Ne ix 
r 
i 
f 

Fexix 
Fexix 
Fe xvn 


12.14 ± 0.0004 
cf. Sect. 4. 2. 1 
13.46 ± 0.003 
13.56 ± 0.011 
13.71 ± 0.004 
13.52 ± 0.003 
13.79 ± 0.007 
13.84 ± 0.07 


0.021 ± 0.0008 

0.022 ± 0.0022 
0.022 ± 0.0022 
0.022 ± 0.0022 
0.022 ± 0.0022 
0.022 ± 0.0022 
0.022 ± 0.0022 


2481.48 ± 56.31 

665.03 ± 37.21 
146.65 ± 31.74 
385.38 ± 30.83 
588.98 ± 39.82 
147.29 ± 28.81 
114.18 ± 27.80 


6510 
6000 


2.63 ± 0.61 
[2.62 ± 0.61] (2) 
[fit constrained t 


0.80 
o G=0.8] 


24.96 

26.16 
26.23 
26.29 
26.21 
26.31 
26.32 


O vin 
O vnr 
i 
f 


18.9701 ± 0.0004 
21.62 ± 0.016 
21.82 ± 0.021 
22.11 ± 0.022 


0.0262 ± 0.0004 
0.022 ± 0.002 
0.021 ± 0.004 
0.020 ± 0.004 


2882.96 ± 57.81 
262.49 ± 22.6 
128.77 ± 18.7 
120.9 ± 18.0 


3150 
2439 


0.94 ± 0.2 
[0.94 ± 0.2] (2) 


0.95 ± 0.16 
[0.97 ± 0.16] (2) 


24.29 
15.58 
15.34 
15.32 


N vii 
N vir 

i 
f 


24.8 ± 0.0012 
28.81 ± 0.007 
29.11 ± 0.011 
29.55 ± 0.011 


0.03 ± 0.0012 
0.037 ± 0.008 
0.058 ± 0.012 
0.021 ± 0.007 


1119.05 ± 38.38 
141.33 ± 21.10 
188.23 ± 25.53 
37.5 ± 14.31 


1700 


0.20 ± 0.08 
[0.21 ± 0.09] (2) 


1.60 ± 0.37 
[1.65 ± 0.38] (2) 


15.23 
13.57 
13.57 
12.76 


Fexxi density diagnostics (cf. Sect. 


6.1 


) 


F(A) W 
F(128.73A) 


ISM 






97.87 

102.22 

117.505 

121.22 

128.73 


0.049 
0.043 
0.058 
0.055 
0.058 


67.20 ± 14.50 
94.11 ± 16.1 
92.10 ± 16.00 

< 36 
266.32 ± 20.90 


334 
350 
225 
117 
73 


0.12 ± 0.03 
0.18 ± 0.07 
0.20 ± 0.04 
<r n no 

1 


0.9206 
0.9111 
0.8734 
u.ouoo 
0.8420 


7.17 
6.64 
6.13 

o.oo 

3.64 


Neix consistency check (cf. Sect. £.jyj 


) 








Fexix 
Fexix 
Fexix 


14.66 ± 0.007 
101.63 ± 0.007 
108.45 ± 0.007 


0.02 ± 0.002 
0.043 ± 0.006 
0.054 ± 0.006 


< 80 
77.67 ± 14.55 
203.53 ± 19.42 


5010 
350 
311 




0.999 
0.9124 


26.99 
6.69 


Fexvn optical thickness (cf. Sect, 


F(15.03A) 
F(15.27A) 


F(15.03A) ( 2 ) 
F(15.27A) 




Fe xvn 
Fe xvii 


15.026 ± 0.001 

15.27 ± 0.008 


0.021 ± 0.001 
0.023 ± 0.003 


1018.44 ± 38.93 
364.71 ± 28.92 


4500 
4500 


2.79 ± 0.25 


[2.81 ± 0.25] (2) 


27.21 
27.42 


Fe temperatures (cf. Sect, pi) 


A[cts/cm 2 ] (2) 






Fe xvii 
Fe XVIII 

Fe xxn 
Fexx 
Fexx 


15.28 ± 0.004 
16.082 ± 0.0031 
117.25 ± 0.093 
118.81 ± 0.014 
121.98 ± 0.009 


0.024 ± 0.003 
0.03 ± 0.004 
0.053 ± 0.06 
0.0623 ± 0.013 
0.0675 ± 0.009 


375.73 ± 29.30 
572.78 ± 35.95 
483.14 ± 25.91 
97.87 ± 16.98 
180.88 ± 19.23 


4500 
4500 
225 
225 
138 


14.22 ± 1.11 
21.07 ± 1.32 
78.56 ± 4.21 
16.26 ± 2.82 
34.19 ± 3.64 




26.42 
27.19 
6.15 
6.02 
5.29 



^' Fluxes accounting for A c g and ISM as in last but one column. ' 2 ' Fluxes corrected for A c e listed in the last column. 



